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The vortex core structure in a d-wave superconductor is analyzed on the 
basis of the quasi-classical Eilenberger theory beyond the Ginzburg-Landau 
framework. The current and magnetic field distributions around an isolated 
vortex break circular symmetry seen in s-wave pairing and show four-fold 
symmetry, reflecting the internal degrees of freedom in d-wave pairing, i.e. 
kl — ky in reciprocal space through the low lying quasi-particle excitations. 
The peculiar orientation of the flux line lattice observed recently in a cuprate 
is argued in light of the present theory. 
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Much attention has been focused on the high temperature cuprate su- 
perconductors. Although precise pairing symmetry has not been identified 
yet, it is recognized that d^.^-wave symmetry is most probable. 

Recently, Keimer eial.[Q reported the small angle neutron scattering 
experiment on YB^CusOy to see the flux line lattice (FLL) in a magnetic 
field region; 0.5T < H < 5T and found that the vortices form an oblique 
lattice with two equal lattice constants and an angle of 73° between the two 
primitive vectors. FLL is oriented such that the nearest neighbor direction of 
vortices makes 45° from the a-axis (or o-axis) of the underlying orthorhombic 
crystal lattice when the field H is applied to the c-axis. This result is highly 
non-trivial to be understood within the framework of the Ginzburg-Landau 
theory (GL). In order to stabilize such an orientation of FLL, one of the 
primitive vectors must be aligned, making an angle of 9° from the a-axis. 
This means that we need the orientational energy of a higher order form 
Aq cosnO with n = 10 in the GL free energy, where 9 is the angle between 
one of the primitive vectors and the a-axis (or o-axis). This is quite unlikely 
in the spirit of the order parameter (A ) expansion. 

Within the framework of GL, which is widely and conveniently used 
to discuss electromagnetic properties of cuprates, s-wave and d-wave super- 
conductivity are not distinguished as far as the order parameter is single 
component because only the number of the order parameter components are 
relevant to construct the GL free energy. This is particularly true for the 
vortex structure. Therefore, the vortex core structure of d-wave pairing has 
been considered to be cylindrically symmetric with a singularity at the cen- 
ter, as in s-wave superconductors, within GL framework. The supercurrent 
and magnetic field distributions around an isolated vortex line are also cylin- 
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drically symmetric unless the crystalline anisotropy is implemented. (Note 
that the Cu-0 plane in cuprates is quite isotropic in this respect.) 

Motivated by the neutron experiment and the above considerations, we 
undertake to examine the vortex core structure of d-wave superconductiv- 
ity. This study underlies the fundamental understanding of electromagnetic 
properties of cuprates. We hope that the internal degrees of freedom of a 
d-wave pair characterized by the relative coordinates such as A(k) (x.k x — k y 
might affect on determining properties related to spatial variations which are 
governed by the center of mass coordinates of a pair, such as a core structure 
or interface 0] where the pair function is forced to vary spatially. 

The Eilenberger theory|| is best suited for our purposes. Because this 
quasi-classical approach was originally designed to describe the conventional 
s-wave superconductors beyond the GL framework microscopically. This 
method easily allows us to treat the two different coordinates systems men- 
tioned above and to discuss the physics of the length scale ~ £o (the coherence 
length), ignoring the shorter scale kp 1 (kp\ Fermi wave number). 

The purpose of this paper is to analyze the vortex structure for a d- 
wave pairing beyond the GL framework, namely, within the quasi-classical 
Eilenberger theory. We will conclude that reflecting the internal degrees of 
freedom of a d-wave pairing function, the vortex core structure is anisotropic 
and the cylindrical symmetry around a vortex line seen in s-wave pairing is 
spontaneously broken. The associated supercurrent distribution and mag- 
netic field distribution also become four-fold symmetric. This spontaneous 
symmetry breaking of the vortex structure is generic in anisotropic supercon- 
ductors. Other non-trivial pairing states which are known to exist in heavy 
Fermion superconductors such as UPt 3 || and superfluid 3 He0 must give 
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rise to own particular anisotropic vortex core structures, corresponding to 
their particular symmetries both of orbital and spin structures of the pair- 
ing function. The cylindrical symmetry around a vortex line is realized only 
when the pairing is s-wave. 

We starts with a transport-like Eilenberger equation for the quasi-classical 

n . \ Y iuj n — A(r,9)\ ^,-> n ■ \ i 
iv F ■ Vg(r,9,iu n ) + [(^ A */ f m -iu n ) > 9{r,9,tuj n )\ = 

supplemented by the normalization condition: g(f, 9, iu n ) -g(r, 9, iu n ) = —n 2 ■ 
1. The Fermi velocity is vf (The Fermi surface is assumed to be a circular). 
The order parameter A(r, 9) is described by two kinds of the coordinates, 
namely, the two-dimensional center of gravity coordinates r = (x, y) and the 
relative coordinate 9 which is measured from the a-axis (or x-axis). The 
bracket [ , ] is a commutator. The self-consistent equation is given 

/jnr 
—V(9,9')f(f,9',iu n ) 

where the pairing interaction V(9, 9') is assumed to be separable. The density 
of states at the Fermi surface is iVo- The associated supercurrent j{r,9) is 
given in term of g(r, 9, iui n ) by 




The magnetic field is determined through —j — V X H. Instead of solving 
the above Eilenberger equation, it is more convenient to use the following 
parametrization devised by Schopohl and Maki|| 
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, _ _ . 1 / 1 - ab 2ia \ 
9 ~~ m \ + ab \ -lib -1 + ab) ■ 

The Riccati equation is satisfied by a(r\\) and b(r\\): 

The other unknown quantity b(r\\) is related by b(r\\) = —a(—r\\)e~ 2ie by 
symmetry in the isolated vortex case under consideration. Here we have 
taken the coordinate system: u = cos 9x + sin 9y,v = — sin 9x + cos 0y, thus 
a point r = xx + yy is denoted as f = r\\u + r±v. 

We take the following form for the d-wave vortex structure as a test 
potential: 

A(r, k, 6) = A tanh^ • e** cos 2^ 

and V(0, 6 1 ') = g a cos 29 cos 29' with g d being the coupling constant for the 
d-wave channel. The phase factor is given by e 1 ^ = 1= . The numerical 

y/x 2 +y 2 

computation has been done using the standard algorithm for solving the first 
order differential equation. 

The stereographic view of the current distribution near the vortex core 
is shown in Fig.l. It is clear that the amplitude of \j(f)\ is not circularly 
symmetric. The breaking of the cylindrical symmetry is stronger in the core 
region while toward the outer region the cylindrical symmetry is recovered 
gradually. The four-fold symmetry of \j(r)\ arises from the d-wave nature of 
the relative coordinate of the pairing function A(k) oc k^ — ky. The four small 
peaks around the core move inside toward the core when the temperature 
(T) is lowered, indicating that the core region is narrowed with decreasing 
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T (compare Figs. 1(a) and (b)). This feature coincides qualitatively with a 
general rule known for conventional s-wave superconductors and confirms 
an interpretation that it arises from Fermionic excitations bounded by the 
core associated with d-wave pairing by Volovik[[ITJ. The current distribution 
of \j(r)\ extents toward the 45° direction from the x-axis and its equivalent 
directions. This is physically because there exist the local quasi-particle exci- 
tations along those directions which contribute to inducing the supercurrent. 
In fact the local density of states of the quasi-particle excitations || has a 
characteristic four-fold symmetric shape for low energies (< A ) and almost 
circular shape for higher energies (> Aq). The latter is responsible for almost 
circular distribution of \j(r) \ in the outer region. It is known that in the core 
region physical quantities are dominated by lower energy excitations ||1 1|. 

The associated magnetic field distribution H(r) is displayed in Figs. 2(a) 
and (b). It is evident from Fig. 2(a) that H(r) is also four- fold symmetric. 
The field extends along the x-axis and its equivalent directions, rotated by 
45° from the current distribution. This is because the field is much screened 
by the induced current along the 45° direction, reflecting that the local quasi- 
particle excitations with lower energies deviates from the circular shape. The 
field distribution becomes almost circular beyond a characteristic length scale 
£o- We show a stereographic view of H(f) in Fig. 2(b) where the ridges are 
barely seen along the and its equivalent directions. 

While our calculation is based on the test potential as a first step, 
the feature of our results reflects qualitatively the essential nature of the 
d-wave pairing. The quantitative aspect might be modified by performing 
the self-consistent calculation. This is particularly true for the degree of the 
anisotropy in H(f). In fact, we can certainly anticipate that the anisotropy 
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in H(f) here is further amplified by a self-consistent computation. Because 
the resulting order parameter profile around the core, calculated by assuming 
a circular symmetric form: A(f) oc tanh| • e 1 ^ as a starting test potential, 
comes out to be a anisotropic one (Note the square shapes near the center in 
the contour plot of Fig. 3). This will further lead to the anisotropic current 
and field distribution. The full self-consistent computation is left for a future 
problem. 

Finally, we point out a possibility that other components of the order pa- 
rameter may be induced around a d-wave vortex: As shown in Fig. 4 as an ex- 
ample, the s-wave component evaluated by A s (f) = N g s 27rT ^2 J 0', iu n ) 
is induced only in the core region of a d-wave superconductor if the pairing 
interaction has a small s-wave component g s . Associated with the spatial 
inhomogeneity of the order parameter at a vortex or interface, other com- 
ponents may be induced. This tendency is in agreement with the similar 
quasi-classical analysis of interfaces [|J and with the GL theory analysis of 
vortices |T^|. The mixing of two kind of order parameters, e.g. s + i<i-pairing 
vortex further breaks spatial symmetry, giving rise to the two-fold symmetric 
local density of states. 

Based on our calculations, let us argue how the flux line lattice (FLL) 
differs from the conventional superconductors: The FLL in a d-wave could be 
different from the regular triangular lattice in a s-wave because the latter is 
led due to the isotropy of the vortex core. The distorted vortex core structure 
in the d-wave must be taken into account when determining the structure of 
the FLL. 

As for the FLL orientation relative to the underlying crystalline lat- 
tice, in the traditional treatment based on the GL theory the higher or- 
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der derivatives, or the so-called non-local corrections are considered to de- 
rive an orientational energy. Here the orientational energy comes from the 
anisotropic current distribution and field distribution. We can speculate 
a possibility that when placing a nearest neighbor vortex at a certain dis- 
tance from the center vortex, it could be advantageous to place it where the 
magnetic field is low to minimize the condensation energy loss, that is, the 
45° direction in our case. This may explain the observation of the neutron 
experiment mentioned earlier. These arguments should be substantiated 
by the detailed calculations beyond the GL theory. Here we only show the 
essence. Together with these theoretical investigations, the experimental ef- 
forts to understand the fundamental electromagnetic properties of the FLL 
associated with a d-wave pairing are certainly needed, for example, the fur- 
ther small angle neutron scattering or STM experiments are desirable to see 
the actual field distribution, or to directly probe the current distribution 
around a vortex core, which could distinguish it from s-wave pairing. 

We thank N. Schopohl, K. Maki, D. Rainer and K. Takanaka for useful 
theoretical information and B. Keimer for discussions on their experiment. 



8 



References 

[1] See for example, papers in the Proceedings of M 2 S Conference in Physica 
C 235-240, 1 (1994). 

[2] B. Keimer, et.al, Phys. Rev. Lett. 73, 3459 (1994). 

[3] K. Takanaka, Prog. Theor. Phys. 50, 365 (1973). K. Takanaka and K. 
Machida, private communication. 

[4] The order parameter mixing is discussed to occur at interfaces in M. 
Matsumoto and H. Shiba, preprint 

[5] G. Eilenberger, Z. Phys. 214, 195 (1968). See for review, J. W. Serene 
and D. Rainer, Phys. Rep. 101, 221 (1983). An analysis of the vortex 
structure of a s-wave superconductor is done in L. Kramer and W. Pesch, 
Z. Phys. 269, 59 (1974). 

[6] See for example, K. Machida et.al., J. Phys. Soc. Jpn. 64, 1067 (1995). 

[7] M. M. Salomaa and G. E. Volovik, Rev. Mod. Phys. 59, 533 (1987). 

[8] N. Schopohl and K. Maki, preprint. K. Maki et.al, Physica B204, 214 
(1995). 

[9] U. Klein, Phys. Rev. B40, 6601 (1989) and B41, 822 (1989). 

[10] G. E. Volovik, JETP Lett. 58, 455 and 469 (1993) and Physica C235- 
240, 2411 (1994). 

[11] F. Gygi and M. Schluter, Phys. Rev. B43, 7609 (1991). 



9 



[12] P.I. Soininen etal, Phys. Rev. B50, 13883 (1994) and A. J. Berlinsky 
et.al, preprint. Also see, Y. Ren et.al, Phys. Rev. Lett. 74, 3680 (1995). 



10 



Figure Captions 

Fig.l: Amplitude of the current density \j(r)\ in the core region at T/T c = 
0.05 (a) and T/T c = 0.5 (b) in the unit of — 2eNoVF^-o (T c is the transition 
temperature). There are four small peaks around the core in the 45° and its 
equivalent directions. These peaks become clear and move toward the core 
upon lowering T. Note the change of the length scale in (a) and (b) which 
is normalized by £o- 

Fig. 2: Contour plot (a) and stereographic view (b) of the magnetic field 
distribution at T/T c = 0.05. The core region is focused in (a). The field 
extends along the x-axis and ?/-axis. 

Fig. 3: Amplitude of |A(f)|, the k^ — ky component, at T/T c = 0.05 calculated 
by using the resulting Green function. It is normalized by Ao- The amplitude 
is suppressed in the x-axis and y-axis directions as is seen from the square 
shapes in the contour plot. 

Fig. 4: Amplitude of \A s (f)\/(g s / g d ), the s-wave component, at T/T c = 0.05 
calculated by using the resulting Green function. The s-wave amplitude is 
induced in the core region. The contour plot shows that the profile has a 
four-leafed clover shape as coinced with the GL prediction |12[|. 
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Fig. 1 (a) 




Fig. 1 (b) 




Fig. 2 (a) 



Fig. 2 (b) 
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